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Bianchi type I cosmological models are studied that contain a stiff fluid with a shear 
viscosity that is a power function of the energy density, such as ( = ae". These models are 
analyzed by describing the cosmological evolutions as the trajectories in the phase plane of 
Hubble functions. The simple and exact equations that determine these flows are obtained 
when n is an integer. In particular, it is proved that there is no Einstein initial singularity 
in the models of < n < 1. Cosmologies are found to begin with zero energy density 
and in the course of evolution the gravitational field will create matter. At the final stage, 
cosmologies are driven to the isotropic Fnedmann universe. It is also pointed out that 
although the anisotropy will always be smoothed out asymptotically, there are solutions 
that simultaneously possess non-positive and non-negative Hubble functions for all time. 
This means that the cosmological dimensional reduction can work even if the matter fluid 
having shear viscosity. These characteristics can also be found in any- dimensional models. 
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1 Introduction 



The investigation of relativistic cosmological models usually has the energy momentum ten- 
sor of matter as that due to a perfect fluid. To consider a more realistic model one may take 
into account the dissipative processes that are caused by the viscosity and that have already 
attracted the attention of many investigators. Misner [1] suggested that the anisotropy in 
an expanding universe would be smoothed out by the strong dissipative process due to the 
neutrino viscosity. Viscosity mechanisms in the cosmology can explain the anomalously high 
entropy of the present universe [2,3]. Bulk viscosity associated with the grand- unified-theory 
phase transition [4] may lead to an inflationary scenario [5,6]. (The inflationary cosmology 
invented by Guth [7] in 1981 is used to overcome several important problems arising in the 
standard big bang cosmology.) 

Murphy [8] obtained an exactly soluble isotropic cosmological model of the zero-curvature 
Priedmann model in the presence of bulk viscosity. The solutions Murphy founded exhibit the 
interesting feature that the big bang type singularity appears in infinite past. Exact solutions 
of the isotropic homogeneous cosmology for the open, closed, and flat universe, were found 
by Santos et al.[9] when the bulk viscosity is the power function of energy density. However, 
in some solutions, the big bang singularity of infinite density occurs at finite past. It is thus 
shown that, contrary to the conclusion of Murphy, the introduction of bulk viscosity cannot 
avoid the initial singularity in general. Anisotropic models with bulk viscosity, which is the 
power function of energy density, have been discussed in detail in our previous papers [10,11]. 

Belinskii and Khalatnikov [12] presented a qualitative analysis about Bianchi type I 
cosmological models under the influence of shear viscosity. They then found a remarkable 
property that near the initial singularity the gravitational fleld creates matter. Recently, 
Banerjee et al.[13] obtained one Bianchi type I solutions for the case of stiff matter by using 
the assumption that shear viscosity (j]) is the power function of the energy density (e), i.e., 
1] = ae", where a is a constant. However, Banerjee et al. merely analyzed the behavior of 
the cosmological models for some values of n. 

In this paper we will investigate the cosmological models again. Since we study these 
models by describing the evolution of cosmologies as the flow in the phase plane of Hubble 
functions, we can clarify the property of the models with any value of n. In particular, 
we prove that there is no Einstein initial singularity in the models with < n < 1. The 
cosmologies have zero energy density in the initial phase and then the shear viscosity causes 
the gravitational field to create matter during the evolution. At the final stage, cosmologies 
are driven to the isotropic Priedmann universe. Although the anisotropy in the universe 
is smoothed out asymptotically, we point out that there are solutions that simultaneously 
possess non-positive and non-negative Hubble functions for all time. In view of this fact. 
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we then consider the higher-dimensional theory and show that the cosmological dimensional 
reduction may work in cosmological models which have shear viscosity. The models extended 
to higher dimensions have also been analyzed according to the methods described in the 
present paper, the results show that they all share the same characters. 

The organization of this paper is as follows. In Sec. II the derivation of two dynamical 
evolution equations of expansion and shear scalars, which was presented in Rcf. [13], is 
summarized for a convenient reading. In Sec. Ill we consider the axially symmetric Bianchi 
type I models in which there are only two Hubble functions. The evolutions of the cosmology 
are described as the flows in the phase space of the Hubble functions. The characteristics 
of the evolutions of the cosmological models for any n arc then clarified. The Bianchi type 

1 models with three Hubble functions are discussed in Sec. IV, where we also consider 
the higher-dimensional models. Section V is devoted to a summary. In the Appendix the 
solutions in the early stage for the n > 1 models, which explicitly show how the energy 
density approaches zero at the initial singularity, are given. 

2 Einsterin Field Equation 

We consider the (1 + D)-dimensional Bianchi type I space-time with 

D 

1=1 

The energy-momentum tensor for a fluid with shear viscosity is 

r^i/ = (e + P)ut,u^ + pg^u - rjK'fiu, (2.2) 
P^p-{2/D)rju^^, (2.3) 

K-fMu = u^.y + u^.^ + u^n^Uy-x + u^n^u^.x, (2.4) 

where e is the energy density, p is the pressure, and 1] is the shear viscosity, respectively. 
Choosing a comoving frame, where u'^ = 6^, the explicit form of the Einstein equations is 

^ij.i' ~ 2^tiyR = T^u, (2.5) 

for the metric (2.1) and the energy-momentum tensor (2.2) can be written as 

[{D -1)/2D]W'^ -a^ = 6, (2.6) 



the line element 

(2.1) 
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= P + ^ - 2r,H,, (2.7) 



where the Hubble functions H,, the expansion scalar W, and the shear scalar are defined 
by 



(2. 



The trace part of Eq. (2.5) leads to 



2— + -H^2^2a2 
at L) 



D-1 



[e-Dp]. 



(2.9) 



Using relation (2.6) to eliminate the in Eq. (2.9) we obtain 



dW 
"df 



D 



D-1 



[e-p). 



(2.10) 



As a consequence of the Bianchi identity we have 

de 



dt 



+ {e+p)W -Arja^ ^0. 



(2.11) 



Equation (2.6) can yield a relation 



d{ayW^) _ d{e/W' 



dt 



dt 



(2.12) 



After substituting the expressions de/dt and dW/dt in Eqs.(2.10) and (2.11) into the rhs of 
Eq. (2.12) we finally obtain the evolutional equation of a shear scalar: 



dja^W^) 
dt 



= -(— ) 



2D 



D-V W 



+ 4?7 



(2.13) 



When the universe is filled with stiff matter, i.e., e = p, we can, from Eq. (2.10), find the 
solution of an expansion scalar: 



W = 1/t. 



(2.14) 



Using relation (2.14), Eq. (2.13) can be written as 
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(2.15) 



Equation (2.15) can be integrated exactly if n is an integer, as has been shown by Banerjee 
et al. [13]. However, one can only analyze the models for some values of n as a result of 
the fact that these integrated relations are too complex. To overcome these difficulties, in 
Sec. Ill we will analyze the cosmological evolutions by describing them as flows in the phase 
plane of Hubble functions. 

3 Solution in the Phase Plane 

We first consider the axially symmetric Bianchi type I model in the 1+ 3 dimension; the 
Hubble functions therein are denoted as 

H, ^h, H2^H^^ H, (3.1) 

in terms of which one can, from Eqs. (2.6) and (2.8), obtain 

W^h + 2H, a'^^Uh-Hf, e^H{2h + H). (3.2) 

We express h and H in terms of the variables r and Q: 

h = rsinO, H = rcosO. (3.3) 

Thus 

W — r{sin6 + 2cos6), 

y^^[{l-tane)y{2 + tane)% (3.4) 



and Eq. (2.15) can be written as 

/ f y ^ -^Oi [ t-^^dt = + C = ^^W^2n-l ^ c 

J y{\-yY J 2n-l 2n-l 

An 

= — — -r2"-i(sm^ + 2cos^)2"-i + C, (3.5) 
2n — 1 

where C is an integration constant. Because y is the function of the variable 9 only, Eq. (3.5) 
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tells us that we can express r as the function of the variable 6. It is through this property 
that one can study the cosmological evolutions by analyzing these trajectories in the phase 
plane. 

A. Evolutions of the Cosmologies 

1. Fixed point 

Wc only consider the plane where the trajectories are in the regions that satisfy the 
dominant energy condition [14], i.e., e > 0. The evolutions of the cosmology are therefore 
confined to the regions H > and H + 2h > 0, as Eq. (2.6) shows. The evolutions of the 
cosmology should start from a fixed point or infinity and then end in another fixed point or 
infinity in the phase plane. From Eq. (2.13) we know that the fixed points should be those 
with zero energy density, which are at H = ot H + 2h = 0. 

Furthermore, from Eq. (2.13) we can obtain 

^ = -2W - 4^. (3.6a) 

Since rj >0 and W = 1/r, we then obtain a relation 

< r^. 3.66) 

Equation (3.6b) tells us that the shear scalar cr^ should be a decreasing function and that 
anisotropy will be smoothed out asymptotically. Therefore, the original point is an attractive 
fixed point. 

2. Invariant Lines 
Equation (2.15) can lead to 

| = (3.7) 

Substituting relation (3.4) into Eq. (3.7) we then obtain 

^ = ^r2"(cose - sine)(2cose + sine)[{cose + 2sine)cose]''. (3.8) 

From the theory of nonlinear differential equations we know that the zeros of Eq. (3.8) give 
invariant lines in the phase plane of /i x if. They are as follows. 

(i) The isotropic state is cos9 — sinO — 0. In practice, this shall only refer to the 
original point, resulting from the fact that there is no shear viscosity in the isotropic state. 
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Furthermore, since the shear scalar is a decreasing function, the original point must be the 
attractive state. 

(ii) The states 2cosd + sind = have negative energy density and are thus neglected. 

(iii) The states cos6 = or cos9 + 2sin6 = have zero energy density. Because the shear 
scalar is a decreasing function, the vacuum states must be the initial states. 

Using results (i)-(iii) and the exact solution of the n = model (note that we can obtain 
the exact solution for 2n = integer models) we therefore conclude that the cosmologies 
should begin with zero energy density in the initial phase; then the shear viscosity causes the 
gravitational field to create matter during the evolution; and at the final stage, cosmologies 
are driven to the isotropic Friedmann universe. 

B. Singularity 

Equation (3.2) tells us that the energy density can become infinity only if H and/or h 
are infinite. Equation (2.9) tells us that dW/dt can become infinity only if H and/or h are 
infinite. Since the Riemann scalar curvature can be expressed as 

R^2^ + W^ + 2H^ + h'', (3.9) 
at 

one then sees that R can become infinity only if H and/or h are infinite. Therefore, the 
Einstein initial singularity can arise only if H and/or h are infinity, i.e., r ^ oo. 

From Eq. (3.5) we know that r can become infinity only if the integrated value on the 
Ihs is infinity. Since in this section we have discussed that only those states with zero energy 
density may be the infinite value of r, we therefore need only to consider the vacuum states. 

When y —■ ^, then 

^ " ' ^ (3.10) 



t/(l/3-y)- J (1/3-2/)" 

The integration (3.10) is infinite only if n > 1, as easily seen. Thus we have shown that 
models with < n < 1 can avoid the Einstein initial singularity. 
C. Examples 

To give some examples we present the following explicit solutions. 
For n = 1/2, 



J3(l + 2tand) - (2 + tane) 

[r{sin9 + 2cos9)]^'^/^ = c \ =^ -, (3.11) 

/3(1 + 2tan9) + {2 + tan9) 



where C is an integration constant. Using Eq. (3.11) we can plot the trajectories in the 
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phase plane; they then determine the cosmological evolutions. One see that r is finite when 
cos9 = or 1 + 2tan9 — 0; thus this model has no initial singularity. 
For n — 1, 



In 



{l-2taney _ 4a 
(1 + 2tan9) ~ T 



r{sin9 + 2cos0) + C, 



(3.12) 



where C is an integration constant. Using Eq. (3.12) we can plot the trajectories in the 
phase plane; they then determine the cosmological evolutions. One see that r is infinite when 
cos9 = or 1 + 2tan9 — 0; thus this model begins with an initial singularity. 

4 More General Solutions 

We now consider the Bianchi type I models with multiple Hubble functions. Prom Eq. (2.7) 
we can obtain 



where C, are the integration constants. Relation (4.2) tells us that one can express all other 
D — 2 Hubble functions in terms of only two Hubble functions. Through the same procedures 
as those described in Sec. Ill we find that the evolutions of the 1 + /^-dimensional Bianchi 
type I cosmological models containing stiff fiuid with shear viscosity as a power function 
of the energy can also be expressed as fiows in the phase plane of H x h... One can also 
show that the characteristics of the cosmological evolutions are just those described in Sec. 
Ill, i.e., cosmologies will begin with zero energy density; in the course of evolution the 
gravitational field will create matter; and finally, cosmologies are driven to the isotropic 
Priedmann universe. In the same way, we can also prove that there is no Einstein initial 
singularity in models with < n < I. 

To give an illustration we consider the (1 + 3)-dimensional model with three Hubble 
functions. We now have the relations 




(4.1) 



Equation (4.1) yields the solutions 



Hi^{l-Ci,)H, + QH2, i 



(4.2) 



//g = (1 - C)Hi + CH2, W^{2- C)Hi + (1 + C)H, 



2, 
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a^^[{l-C + C^)/3]{Hi-H2f. 



(4.3) 



Expressing Hi and H2 in terms of the variables r and 6, 

Hi = rcos9, H2 = rsinO. (4.4) 

Then Eq. (2.15) can lead to 

^ = — r2"(l-C'+C2)(cos^-s«n^)[(2-C)cos^+(l+C')sm^][(l-C)cos2^+2sm^cos^+Csm2^]". 
at 3 , ^ 

(4.5) 

The zeros ofEq. (4.5) give invariant hnes in the phase plane. They are as follows. 

(i) The isotropic state is cos9 — sinO — 0. As discussed in Sec. Ill, this shall only refer to 
the original point, which corresponds to the final state. 

(ii) The states (2 — C)cos9 + (1 + C)sin9 — can be shown to have negative energy 
density and are thus nelected. 

(iii) The states (1 — C)cos^9 + 2cos9sin9 + Csin^9 — are the vacuum states. Because 
shear scalar is a decreasing function, these states must be the initial states. 

Results (i)-(iii) show that the cosmologies will begin with zero energy density and then 
end in the isotropic Fnedmann universe. One can also use the method described in Sec. 
Ill to prove that models with < n < 1 have no Einstein initial singularity. Finally, we 
want to mention that although the amsotropy is smoothed out asymptotically [as Eq. (3.6b) 
shows] and cosmologies will be driven to the isotropic Fnedmann universe eventually, there 
indeed are solutions that possess both nonpositive and non-negative Hubble functions. This 
implies that the cosmological dimensional reduction can work in cosmological models which 
have shear viscosity. This property has been found in our previous paper. Here we show an 
example in order to complete our discussions of cosmological models. 

Consider the n = 1 five-dimensional theory with two Hubble functions of h and H 
corresponding to those of three-space and extra space, respectively. Using Eqs. (2.14) and 
(2.15) we obtain 

W = t-\ a" = ^(l/t')[l + Ce-=^"/2*]-^ (4.6) 
8 

where C is a positive constant. With the definitions of W and cr^ in Eq. (2.8) we then obtain 

h = (l/4t)(l + [1 + Cea;p(-3a/2t)]-'/2), 
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H 



H 



h 



(l/4i)(l - 3[1 + Cexp{-3a/2t)]-^/^), 
(l/4i)(l - [1 + Cexp{-3a/2t)]-^/^), 
(l/4t)(l + 3[1 + Cexp{-3a/2t)]-^/^), 



(4.7) 



(4.8) 



Solution (4.7), in which h is positive while H is negative for all time, can be found when 
C < 8. We have also checked that the simultaneous existence of non-positive and non- 
negative Hubble functions is found in other (including 1-1-3 3) space-time models with other 
values of n. Thus we have clarified the effects of shear viscosity on the characteristics of 
cosmological evolution. 

5 Conclusions 

In this paper we have discussed Bianchi type I cosmological models with a viscous fluid, 
assuming that the shear viscosity is a power function of the energy density, such as ?7 = 
ae" We have presented a detailed study of these models by describing the evolutions of 
cosmologies as the flows the phase plane of Hubble functions. We have clarifled the property 
of these models with any value of n. In particular We have proved that there are no Einstein 
initial singulanties in models with < n < 1. The cosmologies have been found to begin with 
zero energy density; then the shear viscosiy, causes the gravitational fleld to create matter 
during the evolution. At the final stage, cosmologies are driven to the isotropic Friedmann 
universe. We have also pointed out that there are solutions that possess both non-positive 
and non-negative Hubble functions for all time. In view of this fact, we then considered 
five-dimensional theory and gave a solution that explicitly showed that the cosmological 
dimensional reduction can work in cosmological models which have shear viscosity. Models 
extended to other dimensions can also be analyzed according to the procedures described in 
this paper; the results show that they all share the same characteristics. 

The models considered in this paper are only concerned with stiff matter. The same 
models with other matter fields are certainly of interest and remain to be studied. 



Since the behavior of the vanishing energy density at the initial singularity, which is argued 
in Sec. Ill to be a common property belonging to n > 1 models, seems unusual, we now give 
the solution in the early stage to explicitly show how the energy density approaches zero 
asymptotically. 



6 Appendix: Solution in the Early Stage 
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Let e = p and i:) = 3 in Eqs. (2.7) and (2.9). We obtain 



^ + HiW^lriW-2riHi, (Al) 
at 6 



Prom the discussions given in Sec. Ill we know that n > 1 models will begin along the 
invariant lines H — or 2h + H — 0. Thus from Eq. (Al) we find the approximate solution 
as H ^ 0. (The case oi 2h + H ^ could be analyzed in the same way and will give the 
same conclusion) 

To be consistent with Eq. (2.14), W — h + 2H — , one can define 

h^r^ - 5h, {A2a) 
H = 5h/2, (A2b) 

where < Sh < 1. Using the relation of Eq. (3.2), e = H{2h + H), imphes 

e ^ r^Sh. {A3) 



Let 7] = ae"", we then approximate Eq. (Al) as 



The solution for the n = 1 model is 



Sh « Cr^exp{-4:a/3t), {A5a) 
€^Ct~^exp{-Aa/3t), {A5b) 



whereC is an integration constant. Solution (A5) is consistent with the asymptotic form of 
the exact solution found in Eq. (3.9) of Ref. 13. The solution for the n> 1 model is 

5h « [4a{n - l)/3(2n - i)]V(i-n)in/(n-i)^ (^5^) 
e ^ [Aa{n - l)/3(2n - i)] i/(i-n)^i/(n-i)^ (^6^) 



Solutions (A5) and (A6) explicitly show that ngeql models have vanishing energy density 
at the initial singularity. 
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